arXiv:1505.04982vl [physics.atm-clus] 19 May 2015 


On the analysis of photo-electron spectra 


C.-Z. Gao®’^*, P. M. Dinh®’^, P.-G. Reinhard^, E. Suraud®’^ 

°'Universite de Toulouse, UPS; Laboratoire de Physique Theorique (IRSAMC), 

F-31062 Toulouse Cedex, France 

^CNRS; LPT (IRSAMC), Universite de Toulouse, F-31062 Toulouse Cedex, 

France 

^Institut fiir Theoretische Physik, Universitdt Erlangen, D-91058 Erlangen, 

Cermany 


Abstract 

We analyze Photo-Electron Spectra (PES) for a variety of excitation mechanisms 
from a simple mono-frequency laser pulse to involved combination of pulses as used, 
e.g., in attosecond experiments. In the case of simple pulses, the peaks in PES re¬ 
flect the occupied single-particle levels in combination with the given laser frequency. 
This usual, simple rule may badly fail in the case of excitation pulses with mixed 
frequencies and if resonant modes of the system are significantly excited. We thus 
develop an extension of the usual rule to cover all possible excitation scenarios, 
including mixed frequencies in the attosecond regime. We find that the spectral dis¬ 
tributions of dipole, monopole and quadrupole power for the given excitation taken 
together and properly shifted by the single-particle energies provide a pertinent 
picture of the PES in all situations. This leads to the derivation of a generalized 
relation allowing to understand photo-electron yields even in complex experimental 
setups. 
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1 Introduction 


Photo-Electron Spectroscopy (PES) has grown over decades to one of the ma¬ 
jor tools of analysis of the structure and dynamics of atoms, molecules or solids 
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mm- With the increasing availability and versatility of light sources, studies 
of PES are now found in all areas of molecular physics, in atoms and simple 
molecules |1] as well as in more complex systems such as clusters [5] or organic 
molecules [6j. The archetypal case is that a system is subject to a laser held 
with moderate intensity and sufficiently high frequency, such that all valence 
electron states can be ionized in one stroke. In this one-photon regime, the 
PES delivers a printout of the sequence of single particle (s.p.) energies from 
the occupied states [7115] . But PES do also allow to access dynamical features. 
A typical example is the multiphoton regime in which ionization proceeds via 
absorption of several photons which is often achieved by IR photons at mod¬ 
erate intensity. The PES then exhibits the pattern of Multiphoton Ionization 
(MPI) piQllTO] showing successive copies of the s.p. spectrum. Further increas¬ 
ing laser intensity leads to increasingly complex patterns, i.e., above threshold 
ionization (ATI) [TTIIT^ and strong-held ionization [T3]ll4lll5j . 

Simple laser setups are characterized by one well dehned photon frequency 
Wias- The energies Ej j, of the PES peaks from such sources follow a well known 
rule relating them to s.p. energies Si |9] as 

T ) (1) 

where i stands for the s.p. state from which the electron is emitted and u is the 
number of absorbed photons for ionization (z/ > 1 for MPI). In other words, 
the PES yield y (Ekm) is then interpreted as a sum of Dirac distributions : 

3^(Ekin) O . (2) 

i 

For high laser intensity, relation ([^ has to be complemented by the energy of 
the ponderomotive potential Up which represents the average kinetic energy 
of a free electron in the laser held unns]. For the sake of simplicity, we shall 
conhne our discussion in the following to cases in which Up remains negligible. 

With the advent of a large variety of coherent light sources [MIHIIIS] , it has 
become possible to access more complex dynamical scenarios, involving sev¬ 
eral frequencies and/or pulse combinations, for example, in pump-and-probe 
setups m- The latter ones even provide a time-resolved access to dynam¬ 
ics, at ionic pace with femtosecond (fs) [2111^ down to electronic pace with 
attosecond (as) pulses [23|I24II25P26|I27] . The PES of course remains a highly 
valuable tool of investigation of the dynamics, as it basically provides insight 
into electronic dynamics via ionization characteristics. But the simple rule 
Eq. ([^ takes a more complicated expression in the case of a complex light 
pulse. For instance, in a pump-and-probe experiment where there are possi¬ 
bly two laser colors, one should add to the s.p. energies linear combinations 
of the two laser frequencies to obtain the positions of the PES peaks. This 
is all the more involved when it is in the multi-photon regime for which the 
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number of combinations grows fast. In this case, the identihcation of the PES 
peaks to a certain combination of laser frequencies can become tricky. It is 
thus important to understand how the PES builds up in the course of irradia¬ 
tion processes involving complicated laser pulses. The question is in fact even 
more general. There also exist experiments measuring kinetic energies of the 
electrons emitted after irradiation by a charged projectile [2811^ and a hrst 
theoretical exploration was proposed in [30]. This also provides a PES but now 
with a “photon” of mixed frequencies, as the effect of a charged projectile is 
basically to deliver a short electromagnetic pulse which covers a broad band of 
frequencies. Indeed, the shorter the pulse, the broader the band of accessible 
frequencies. For the sake of simplicity, we continue to call PES such kinetic 
electron spectra from ultra-short pulses, even if the notion of a “photon” is 
rather untypical here. The analysis of the PES obtained in such fast collisions 
can also reveal structures which can be attributed to eigenfrequencies of the 
system, e.g., the Mie plasmon frequency in the case of metal clusters PI32] 
or Ceo [33] . It should also be noted that the intense plasmon of metal clusters 
also shows up in laser-driven PES and can be spotted when scanning laser 
frequency around the plasmon frequency, as was outlined some years ago [ 31 ] . 

All in all, the PES thus reflects the typical frequencies present in the system, 
either those delivered from outside by a laser held or intrinsic ones as mostly 
visible in the case of a fast collision. It is the goal of this paper to investi¬ 
gate in detail this interplay of frequencies. We shall consider various types of 
excitations in order to scan a variety of dynamical scenarios. At the side of 
test systems, we perform the quantitative analysis mostly for a simple case, 
that is the He atom, to maintain clear signatures. The paper is organized as 
follows. Section briehy presents the theoretical framework. Then we analyze 
a few typical cases involving various laser setups including state-of-the-art at- 
totrains in Section]^ Once identihed the limitations of the simple rule Eq. Q, 
we propose an alternative, more general, rule to cover any dynamical scenario 
in Section [^ Conclusions and perspectives are hnally drawn in Section 


2 A brief sketch of theory 


2.1 Basics 


Our computations are done with Density Functional Theory (DFT). We work 
here in real-time Time-Dependent DFT at the level of the Kohn-Sham (KS) 
picture [ 33 ] • The system is then described by a set of single particle wave func¬ 
tions which follow the time-dependent Kohn-Sham (TDKS) equations 
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.dipi{r,t) 
' dt 


-v' + v;ff(r) 




(3) 


We use Rydberg atomic units {h = e = 2me = 1) throughout the paper except 
where specihcally stated differently. The term Ws in Eq. (|^ is the KS effective 
potential. It is composed of three parts We = idon + Rh + i4c[^]- The term 
Vion is the ionic background potential which is described by a pseudopotential, 
a simple local one for sodium [37] and a Goedecker-type one for helium [38] . 
The term Vn is the Hartree contribution and the last term I4c represents the 
exchange correlation potential. The latter term is a functional of the local 
density g{r,t) = and we work here within the time-dependent 

version of the Local Density Approximation (TDLDA) with the parametriza- 
tion of Perdew and Wang [39]. The TDLDA is complemented by an efficient 
self-interaction correction term BDl which allows us to properly describe the 
ionization threshold PI12] . 


2.2 Excitation mechanisms 


In this paper, we will use various irradiation processes for the excitation : 
a fast colliding charged projectile, a single laser pulse, a sum of attosecond 
pulses, and finally an involved superposition of IR femtosecond pulse with UV 
attosecond pulses. 

The first type of excitation is an instantaneous dipole boost of all electronic 
wave functions at f = 0, mimicking the collision with a swift ion, that is, 
</9j(r, f = 0) = where D is the dipole operator, t] is the boost 

momentum and is the ground state Kohn-Sham wave function for state i. 

The single laser pulse is modeled as 


Eias(t) = Eo frit) sin((:niast) , (4a) 

f^(t) = cos^ ~ 

where 6 is the Heaviside function, T the pulse duration and cuias its frequency. 
T takes typical values of some hundreds of fs, while we will consider in this 
paper either an IR frequency (cuias = cjir = 0.115 Ry) or a UV one (cjias = 
cuuv = 1-5 Ry in Figure]^ and 0.5 Ry in Figure]^. 

As an example of a rather complex pulse, we will consider an attosecond train 
of UV pulses, that is, a sum of pulses similar to that defined in Eqs. @ but 
with cjias = i^atto in the UV range, a pulse duration T of a fraction of fs, and 
each UV attopulse shifted in time. In other words, the attotrain held reads : 
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^atto(^) = eJ*"**"^ 9{t)fT{t - to) sin(a;atto^) , 
a=0 

ta -|- Ot Ttrain 5 

,(;)-cxpf 

sm exp(^ 0^/(4 In 2) ) ' 

T=iVr,„„/2. 


(5a) 

(5b) 

(5c) 

(5d) 


Some additional parameters have been introduced here, let us explain them. 


In most experiments 


an IR fs pulse is used to generate high 


harmonics which then serve to produce the coherent attotrain, which is usu¬ 
ally modulated by an envelop that we denoted by g in Eq. (5c). Here, we take 

(But other 


in our simulations a Gaussian envelop as in the experiment of 
theoretical calculations use a cos^ envelop instead for a comparison with an¬ 
other experiment |13]-) One can also play on the delay between the IR pulse 
and the attotrain. This delay Ar is thus entering the dehnition of the envelop 
g but also in the shifting times ta in Eq. (5b). In the present work, we £x it 
at Ar = 46 fs. As for Ttrain, h corresponds to the time separation between the 
maxima of two successive attopulses. Finally, the number of UV pulses in the 
attotrain is denoted by N and is typically about 10. It is obvious that Eq. (5a) 
represents a considerably complicated pulse and that it will provide a critical 
test case for our analysis. Figure [T] summarizes the various pulses we will use, 
as it displays the time prohle of a single IR fs pulse (top panel), see Eqs. (|^, 
and that of an attotrain (bottom panel), see Eqs. ([^. 


Finally, we will consider the superposition of an IR fs pulse and a UV attotrain, 
as encountered in experiments [2^27ll43ll44ll45j . This case will represent the 
most difficult excitation mechanism to be understood. 


2.3 Numerical treatment 


The TDLDA equations are implemented on a grid in coordinate space. In the 
present (principle) study, we consider spherical systems so that we can recur 
to a 2D cylindrical representation. This considerably reduces the computa¬ 
tional expense. The ground state of the system is determined by the damped 
gradient method. The time propagation uses T-V splitting. To describe ion¬ 
ization, we use absorbing boundary conditions. For details of the numerics, 
see iHH] . A warning is in order here. Experimental PES and related effects 
belong to electrons in the continuum. It is known that a numerical description 
of outgoing waves requires careful choice of grid parameters to avoid artifacts 
from unwanted discretization of the continuum |18] . We therefore choose here 
very large numerical boxes to avoid such artifacts : for the metal cluster Nag^, 
we use 138.4 ag in the longitudinal direction (=symmetry axis) and 71.2 ag in 
the radial one (with the same mesh size of 0.8 ag in both directions), and for 
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Fig. 1. Time profiles of typical laser fields considered in this paper. Top : IR femtosec¬ 
ond pulse as defined in Eqs. @ , with duration T = 105 fs, frequency aiias = 0.115 Ry, 
intensity I = 10^^ W/cm^. Bottom : train of UV attosecond pulses (attotrain), as 
dehned in Eqs. ([^, with Watto = 1-69 Ry, I = 10^^ W/cm^, attopulse duration 
T = 1 fs, delay At = 46 fs, and attopulse separation Ttrain = 1-33 fs. 

the He atom, we use 105.6 ao and 52.8 ao respectively (the mesh size being in 
that case 0.6 ao). 


2-4 Observables 


The spectral distribution of multipole strength is computed from TDLDA in 

the time domain using spectral analysis following the prescription of [^9ll50] . 

For the example of dipole strength, this proceeds as follows. Whatever the 

excitation mechanism (boost or finite-width laser pulse), we record the dipole 

moment D{t) = f d^rz p(r,t) as it evolves over the TDLDA simulation. Note 

that we consider here the dipole along laser polarization axis, denoted by z, 

or, generally speaking, along the excitation direction and/or symmetry axis 

of the system. The time signal 0(t) is Fourier transformed into frequency 

~ 2 

domain, yielding D{u). The dipole power spectrum is then Pd{(^) cc D{u) . 
We consider here the power spectrum rather than the usual dipole strength 
(which is the imaginary part of D{u) [50]) since we compare it to the PES, 
which is also a power spectrum. 


The central ingredient of the paper is the PES, the spectrum of asymptotic 
kinetic energies of emitted electrons. Its computation has been worked out 
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in detail in several papers [M1I51II52II53] ; for a recent extension to the case 
of strong fields, see |H1]. We briefly present the procedure for the sake of 
completeness. It is directly based on the KS s. p. wave functions |31]- We choose 
a “measuring point” tm located far away from the system and just before the 
absorbing boundaries. We record all s. p. wave functions ipi{rM,t) at that point 
all along the simulation time. Because of the large distance from the center of 
the system, we can neglect the KS held and assume a free particle dynamics 
(strong helds, not encountered here, require to consider electron motion in the 
ponderomotive held of the laser [M])- In addition, because we are close to the 
absorbing boundary and far away from the source, only outgoing waves with 
momentum k = kYM/^M will pass the point ym- This allows us to establish 
a revertible relation uj k between momentum and energy. The PES yield 
y^Lrj^ (-I^kin) can then be obtained from Fourier transformation from time t to 
frequency oj of the KS orbitals ipj: 


N 

ynrj^ (^kin) oc ^ I (Pj (r^, Ekin) I ^ (6) 

i=i 

The quantity fIrjM represents here the solid angle related to the direction of 
ym, Ekin) is the time-frequency Fourier transform of ipj{YM,t) and the 

kinetic energy reads Fkm = k'^/2 = u [5l]. Performing such an analysis at a 
dense mesh of measuring points allows us to compute the fully energy- and 
angular-resolved PES piss] . However, since we are interested here purely 
in spectral features, we shall restrict the analysis to a point along the laser 
polarization (or symmetry) axis. 


3 Typical PES interpretation and limitations thereof 


3.1 A simple example 


We start with illustrating Eq. ([^ through a simple test case, namely the Nag"*" 
cluster irradiated by a single UV laser pulse, see Eq. (|^, of frequency wias = 1.5 
Ry, intensity / = 10^^ W/cm^ and duration T = 232 fs. The ponderomotive 
potential is Up= 2.2 x 10“^ Ry, thus negligible. The resulting total ionization is 
iVesc = 0.038, which means that the test case stays safely in the perturbative 
one-photon regime. The s.p. spectrum of the Nag’*' cluster groups into two 
shells of nearly degenerated levels denoted in harmonic oscillator labeling as 
Is and Ip [32]. Figure displays the obtained PES (black line) for this laser 
pulse. The PES shows the typical pattern of one-photon emission. One can 
nicely identify the two peaks corresponding to the two shells. Is and Ip, of 
occupied s.p. levels (mind that, in 2D, the Ip levels are degenerate), shifted 
by the laser frequency, as expected from Eq. ([^. 
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Fig. 2. Excitation of Nag^ by a UV laser pulse with frequency toias = 1-5 Ry, 
intensity of 10^^ W/cm^, and total duration 232 fs. PES (black curve, lower x-axis 
in each panel) is compared with the spectral distribution of dipole power Pn{Lo), 
one shifted by (bottom panel) and the other shifted by eip (top panel) according 
to Eq. (Q. PES is drawn versus the kinetic energy Ekin of the emitted electron and 
exists, of course, only for positive Ekin, while the dipole power Pd shifted by £i 
also extends to negative energies. The vertical dashed line indicates the emission 
threshold at Ekin = 0. PES and Pd are drawn in arbitrary units and are scaled such 
that they have comparable peak heights. 

Figure]^ also shows the dipole power distribution Pd{(^) computed from the 
response D{t) to the given laser pulse. The power spectrum is compared with 
the PES twice, once shifted by the s.p. energy = —0.64 Ry (bottom panel) 
and once shifted by eip = —0.55 Ry (top panel). These Pd{^^ = -Ekin — £i) 
represent the generalization of Eq. (|^ from one dominant laser frequency wias 
to a full spectrum co. Let us now analyze in more detail Pd{oj) (we remind 
that the scale for u appears as the upper horizontal scale in each panel). It 
shows two prominent peaks : the leading, expected, one at wias arising from 
the laser field and a small secondary peak at lower frequency cuMie = 0.2 Ry 
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corresponding to the mnch celebrated Mie plasmon resonance WMie- As it is 
visible in the right parts of each panel of Figure]^ the laser peak of Pd shifted 
by the two Si perfectly coincides with the two PES peaks. One may be snrprised 
by the appearance of a Mie resonance peak because the mismatch between 
i^Mie = 0.2 Ry and wias = 1-5 Ry is huge. Remind however that the laser 
pulse has a hnite width. Therefore, even at cuMie, there remains a faint piece 
of spectral strength in the laser signal which together with the overwhelming 
responsivity of the Mie plasmon produces this sidepeak. Although visible, it 
is naturally suppressed by orders of magnitude and has no impact on the PES 
here. 

Thus far, the picture is straightforward for such a simple laser pulse. The 
analysis will nevertheless quickly grow in complexity when proceeding to more 
involved pulses with a richer spectral pattern, because it has to be performed 
with the two occupied electron shells in the case of Nag^. To simplify the 
picture at the side of the test system, we will concentrate, from now on, on 
the simpler case of an helium atom with one single occupied level (with two 
electrons of opposite spins). It should also be noted that the Helium case is 
especially interesting in our framework as it involves only one (doubly occu¬ 
pied) wave function. For then the self interaction correction which we use [ID] 
exactly removes the self-interaction error at all orders. 


3.2 Dipole boost excitation 


We start the studies on the He atom by hrst considering a simple boost exci¬ 
tation (simulating a collision with a fast charged projectile). Figure shows 
the emerging PES (black) together with the dipole power spectrum Po (red). 
In that case, the simple rule Eq.([^ is by construction not applicable because 
there is no laser around imprinting its frequency cjias onto the process. On the 
other hand, and following the analysis of Section 23 one can still compare 
the distribution of PES and the dipole spectrum PD(a; = Ekin—£is), as is done 
in Figure]^ The sequence of peaks in PES nicely matches a similar sequence 
in Pd when shifting the scales by the s.p. energy of the single occupied state 
Is (which is here identical to the IP). The agreement between both curves 
is striking, primarily concerning the location of peaks and, to a large extent, 
even at the side of the relative heights of peaks. This conhrms the Endings 
of Figure but now in a case where the simple rule Eq. ([^ is by definition 
meaningless. The generalization of the rule ([^ which we read off from Figure 
2] is that the whole distribution of PES yield is strongly related to the, prop¬ 
erly shifted, dipole power, i.e., 3^(Ekin) ^ PDip> = E\an—^is)) both to be read 
as a function of Ekin. This is the conjecture which we now want to scrutinize 
further. 
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Fig. 3. He atom excited by instantaneous dipole boost of 0.001/ao. Photo-Electron 
Spectra (black line and lower horizontal scale on Ekin) and dipole power Pd shifted 
by the energy of the s.p. level, eis = 1.8 Ry, (red line and upper horizontal scale on 
u). The key for Pd indicates that it can be read in two ways, either vs. oj (upper 
x-axis) as usual for a spectrum, or vs. Ekin (lower x-axis) for comparison with PES. 
The vertical solid line denotes the ionization potential (IP). 

3.3 An attotrain excitation 


The case of a simple boost produces a relatively simple dipole spectrum Pd 
with only a few dominant peaks and minor ones (which, nevertheless, show 
up in the PES). It is thus also interesting to see how the above conjecture 
performs in a complicated case involving several laser frequencies on top of 
the eigenfrequencies of the system. This situation is illustrated in Figure 
again in the helium atom. The latter one is irradiated by the attosecond laser 
train dehned in Eqs. ([^. It consists out of iV = 7 individual UV pulses, each 
one of duration T = 1 fs. The time interval between the end of an attopulse 
and the onset of the next one is 0.33 fs. Therefore, we have Ttrain = 1-33 fs, 
corresponding to a frequency entrain = 27 r/(l. 33/0.0484) = 0.23 Ry to the 
system. The lower panel of Figure shows the result for PES and Pd{oj). 
For the sake of completeness, we also show the power spectrum of the electric 
held Patto(i:^) dehned as Patto(i:^) = |Eatto(i:^)P where Eatto(i:^) is the Fourier 
transform of the laser held of the attopulse train i7atto(t)- The dipole power 
spectrum Pd now takes a seemingly complicated “comb” structure which, 
however, can be understood in simple terms. One identihes a wide peak located 
around the attosecond frequency (Uatto = 1-69 Ry (see upper horizontal scale; 
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Fig. 4. Lower: PES (black) and dipole power Pd (red) for a He atom excited by the 
attosecond pulse train (APT) as given in Eq. ([^, see text for details. The green line 
shows the power spectrum Patto of the laser field delivered by the attopulse train. 
The vertical dashed line indicated the ionization threshold. Upper: The optical 
response of the He atom. 


mind that there is a Gaussian envelop to the attotrain, see Eq. (5a), so that 
the maximum of Patto is not exactly at cuatto)- If has a large width of about 0.3 
Ry related to the short duration of only 1 fs. Because the signal is repeated in 
time (attotrain), this broad structure is overlayed by more intense and sharper 
structures precisely separated from each other by cjtrain, corresponding to the 
repetition rate of the pulses. As the IP of the He atom is 1.80 Ry, only the 
high energy wing of these structures shows up in the PES. And the agreement 
between Pd (a;) and PES is striking for this part of the spectrum : Again, both 
the position and relative heights of major PES peaks just map the dipole 
spectrum (mind the logarithmic ordinate scale and the orders of magnitude 
suppression of higher energy PES peaks). 
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It is also interesting to note that the dipole power Pd{oj) (red curve) does not 
follow exactly the power of the external held Patto(i^) (green curve). To un¬ 
derstand these discrepancies, the upper panel shows the pure optical response 
P])°°®t(a;) obtained from a pure dipole boost, as that shown in Figure]^ This 
rehects the dipole eigenmodes of the system. There are marked differences 
between Patto and Pd exactly at the places where we observe strong peaks in 
pboost_ 'jpjg indicates that the spectrum of the dipole response Pd{^) com¬ 
bines structures from both the laser held (i.e. from Patto) and the strongest 
eigenmodes of He in P^°°®*. This by itself is not surprising. What is more in¬ 
teresting is the fact that the PES precisely matches the combination of both 
spectra as given in Pd, down to details of peak splittings. This is especially 
clear for very low energy PES peaks in which secondary peaks of pboost(^) 
still visible. This example of an attotrain thus conhrms our previous hndings, 
namely that Pd(w = Pkin —£^is) is strongly related to the PES. The simple 
estimate remains robust even for an involved pulse plus admixtures of system 
frequencies in Pd- 


4 Towards a general and robust interpretation of PES 


4-1 The surprising case of an attotrain combined to an IR laser 


There remains an interesting case to be explored in connection to attosecond 


trains. In most experiments 


j , an IR pulse is used to generate high 


harmonics which then serve to produce the coherent attotrain. Therefore, the 
attotrain is shot on top of the IR pulse, the latter one provoking slow and 
gentle dipole oscillations of the electron cloud. We use here an IR pulse with 
the prohle Eq.(4a) and the parameters (Uias = i^ir = 0.115 Ry, an overall pulse 
length T = Tir = 105 fs, and a held strength of Eq = 0.0034 Ry/ao (that is, an 
intensity Rr = 10^^ W/cm^). The IR pulse alone does not ionize the ground 
state of He atom, while the UV attotrain on top of the IR pulse leads to a 
measurable ionization. The result of such a setup is shown in Figure]^ again 
comparing the PES and the dipole power Pd- At variance with all previous 
cases, the two spectra do not coincide anymore. Only half of the peaks of the 
PES can be found in the dipole spectrum Pd- Another sizeable ingredient is 
thus clearly missing in this case. 


4-2 Beyond dipole response 


Let us remind that Pd is the Fourier transform of the time evolution of the 
electronic dipole. A laser pulse, being oc D, therefore excites predominantly 
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Fig. 5. PES (black) and dipole power Pjy (red) for a He atom excited by a combi¬ 
nation of an attotrain ([^ (as used in Figure]^ with an IR pulse, see text for details. 
The vertical solid line indicates the ionization threshold. 

the dipole. However, strong pulses can also induce higher order deformation 
of the electron cloud. For example, a dipole shift of the many-electron wave 
function is realized by the operation 

1$,) = e-"^|<Fo) ^ (l - idP - |<ho) (7) 

where P is the operator of total electron momentum and d the size of the shift. 
Only the term oc is usually accounted for in case of very weak helds as 
they are used typically for nanosecond pulses. Such extremely weak pulses can 
excite only modes with odd parity. Pulses in the femtosecond regime, as used 
here, deal with stronger (although still moderate) intensities. And this brings 
also the term oc P^ into play which, in turn, triggers modes with even parity. 
The leading even parity modes are found in the quadrupole (deformation) and 
monopole (stretching/compression) channels. The corresponding observables 
are : 


Q = J dvQQ = J dr (3z^ - - 2 /^)p(r) (8) 

M = y drMp(r) = j drr^p(r) ( 9 ) 

where p(r) is the local electron density. Both can be constructed by angular- 
momentum reduction of a tensor built from the square of a dipole r ® r [HB] • 
Let us consider for the moment the quadrupole case. Considering the deformed 
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state defined in Eq.([^, we find 


($,| Q |$d) = ($o| Q |$o) -id (<J>o| [Q, P] |$o) +y ($o| [A [Q, P]] |<^>o) (10) 


=0 


=0 


7^0 


which shows that the second order term can produce, indeed, a hnite quadrupole 
momentum. Whereas the observable Po can measure only the odd-parity 
modes, the PES measuring operator is an outgoing wave, approximately (^k ot 
which contains odd and even parities. Thus the PES is, in principle, able 
to record also even-parity modes. The idea is now that the PES peaks which 
are missing in Pd correspond to even-parity excitations and thus should show 
up in the power spectrum of the quadrupole Pq and/or monopole Pm- 


4-3 A pure quadrupole excitation as a proof of principle 

Before attempting to analyze the spectrum of Figure|^with even-parity strengths, 
one should check whether such a picture makes sense in a “cleaner” case 
avoiding that the quadrupole contribution has to be hgured out from a mix 
with a dominating dipole. It would thus be interesting to consider a pure 
quadrupole case excluding any dipole. Although it is not clear how to excite a 
pure quadrupole in practice, it provides an interesting Gedanken experiment 
worth being explored. We thus consider an instantaneous quadrupole boost to 
the electron cloud of the He atom, as 

4=,(r) ^ e‘'«4=,(r) (11) 

where A provides the amplitude of the boost. The He ground state has even 
parity and so does the excitation. Thus the dipole moment remains zero 
(within the numerical accuracy) all over the time evolution and no dipole 
eigenfrequency is excited in the process. Note that the quadrupole boost 
mostly leads to quadrupole oscillations, but the pattern are accompanied by 
small monopole oscillations which, however, are marginal for the present case. 

The PES associated to such a quadrupole boost in a He atom is shown in Fig¬ 
ure]^ We also plot on the same hgure the quadrupole power Pq which exhibits 
the actually excited eigenmodes of the system. The coincidence between PES 
and Pq is again remarkable, both in peak positions and relative amplitudes. 
This result tells us that the PES displays the actually excited eigenfrequen- 
cies of the system, whatever way the system is excited. In the usual case of 
dipole excitations, the dipole response is (not surpassingly) exemplihed. In 
a pure quadrupole (no dipole) excitation, the system mostly responds in the 
quadrupole channel and if one performs a monopole excitation, the system 
dominantly responds in the monopole channel (not shown here). The actual 
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Fig. 6. PES and quadrupole power spectrum Pq for excitation of the He atom by a 
quadrupole boost, see Eq. ( pH| , of amplitude of 0.001/ag. 

share of each mode in the response strongly depends on the excitation but 
as well on the system. As a consequence, one can expect that, for a general 
excitation, all channels may be excited (though with different strengths) and 
therefore, one should consider all eigenfrequencies on the same footing. In sim¬ 
ple cases, one multipole is dominant so that the other ones do not show up, 
but this holds not true in general. 


4-4 Back to the combined IR + attotrain case 


To validate the concept for a more involved case, we come back to the He atom 
excited by an IR-I-APT where the dipole response alone could not explain the 
PES, see Figure]^ We now add the quadrupole spectrum Pq to the analysis 
(the monopole spectrum Pm being two orders of magnitude suppressed com¬ 
pared with Pq, it is not considered here). The result is shown in Figure 
First note that both P^ and Pq take similar values after the same rescaling 
to make them match with the PES scale. This is a clear indication that, in 
this case, both dipole and quadrupole channel do respond with comparable 
amplitudes. When piling up both Po and Pq, one hnds a perfect reproduction 
of the PES, the quadrupole spectrum Pq precisely bringing the peaks missed 
by the dipole spectrum. Again the two spectra reproduce both positions and, 
to a large extent, relative amplitudes, of the PES peaks. 
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Fig. 7. The same as Fig. but including quadrupole power spectra Pq. The Pd 
and Pq are scaled with the same factor to match the height of the PES peaks. 

4-5 The Multi-Photon Ionization regime revisited 


Another typical scenario in laser driven dynamics is multi-photon ionization 
(MPI). We investigate this regime here using a simple laser pulse Eq. (4a) 
with frequency of cuias = 0.5 Ry, intensity I = 10^^ W/cm^, and duration 
^puise = 105 fs. The frequency, although in the UV regime, is far below the 
IP of 1.8 Ry. It requires at least four photons to lift an electron into the 
continuum. This represents a typical MPI case. The result is shown in Figure 
It is obvious that the dipole power Pd alone cannot explain the whole PES. 
Remind that a sufficient pulse intensity is required in MPI to allow for multi¬ 
photon processes which, in turn, makes co-excitation of even parity modes very 
likely. The quadrupole power Pq is sizable, but Pd and Pq together cannot 
explain all pattern of the PES in this case. Mind that they both have been 
scaled the same amount for comparison with the PES, but their relative values 
have not been touched (the same holds for the monopole below). In fact, Pq 
is of minor importance here. There is a sequence of PES peaks at low 
which stems clearly from the monopole modes seen in Pm- Thus monopole 
and quadrupole which stay formally at the same level of importance should 
always be included together in the analysis. It is concluded that the pattern 
of PES can be reproduced by a proper combination of monopole, dipole, and 
quadrupole strength distributions. 


To conclude this discussion of the MPI regime let us mention the associated 
regime of above threshold ionization (ATI) induced by high infrared intensities 
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Fig. 8. PES (black), powers spectra of dipole Pd (red), quadrupole Pq (green), 
and monopole Pm (blue) for a He atom excited by a UV pulse, see Eq. Q, with a 
frequency of tcias = 0.5 Ry, intensity 10^^ W/cm^, and duration T = 105 fs. 

with I > 10^^ W/cm^ [TU1I42II57IIS5] . In this case, the ponderomotive energy 
Up induces an energetic shift which can become comparable to the laser fre¬ 
quency : for wiR = 0.11 Ry, we have Up = 0.057 Ry, which is not negligible 
anymore. For the sake of simplicity we avoided, in all the above discussions, to 
consider cases where the ponderomotive shift was important. But in the ATI 
regime we have to account for this effect. This practically does not change our 
analysis but for the fact that the multipole spectral distributions have to be 
shifted hj e +Up instead of e only. But once this precaution has been taken we 
hnd that the present spectral analysis of the ATI PES perfectly holds, which 
perfectly supports our above conclusions. 


5 Conclusions and perspectives 


We have discussed in this paper the relations between a photo-electron spec¬ 
trum (PES) and spectral strength distributions of the basic multipole opera¬ 
tors (dipole, monopole, quadrupole) for irradiated atoms or clusters. In this 
investigation, we used a He atom as a test case to have a simple single par¬ 
ticle (s.p.) spectrum at the side of the system. For the laser pulses, on the 
other hand, we consider a variety of scenarios from a simple boost over a 
one-frequency pulse to complicated mix of IR and atto-second pulses. 
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Traditionally, there is a simple rule of thumb which relates the position of 
peaks in the PES to the energies of the occupied s.p. states plus an appropriate 
multiple of the laser frequency. The simple rule Q relates peak positions 
with s.p. energies, or in other words, the consequent estimate of the PES 
yield, see Eq. ([^, relates peaks with peaks. This rule cannot be applied to 
more complex light pulses which embrace a couple of different frequencies 
as, e.g., in case of fast collisions with a charge ion or of the mixed pulses 
typically used in attosecond physics. A natural generalization of it is to model 
the spectral distribution of PES by the spectral distribution of dipole power 
Pd{oj) augmented with s.p. energies. This approach is indeed able to provide a 
correct picture of PES in some cases. But part of the PES maxima is missing 
in others. A closer analysis shows that the PES is sensitive to excitations with 
odd partity as well as those with even parity. However, the dipole spectrum Pd 
sees only the odd-parity modes. Therefore one must also extend the estimate 
by strengths of even-parity multipoles. The monopole and quadrupole operator 
are the most important in that regime. The idea is thus to include, in addition 
to Pd, the spectral distribution of monopole power Pm and of quadrupole 
power Pq into the generalized estimate, which amounts to 

y{Ekin} Tj M PM {,Pkm T T^DPoiEkin T 'QqPq{^E]^yyi (12) 


where the sum runs over all occupied s.p. states i, and the 77 ’s factors take care 
of the relative strength of each eigenmode. Experience gained in the present 
study suggests that the scaling factors can be taken as being the same for all 
modes, i.e. rjM = f]D = Vq- This spectrally extended rule (12) has been con- 
hrmed for all the types of excitation we have explored, with the contribution 
of the three terms in the rule varing from case to case. This mapping between 
response(s) and PES we thus found turns out to be extremely robust. 


All in all, the emerging picture remains simple in the sense that the PES still 
provides a direct printout of the frequencies contained in the response of the 
system. The simple, traditional, rule ([^ is recovered when exciting by a sim¬ 
ple mono-frequency pulse sufficiently far from any resonance of the irradiated 
system. This is the case where we recognize only the laser frequency (and s.p. 
energies) in the PES. The new rule (12) covers much more general situations. 
A typical example is attosecond physics where a complex mix of IR pulse 
with trains of attosecond pulses is often used. There remain widely disputed 
questions on the mechanisms underlying the response of the system to such 
complex excitations. A typical observable here is the total ionisation which is 
found to oscillate with the delay between IR pulse and attotrain P^27p59j . Our 
preliminary computations show that these modulations are (not surprisingly) 
directly reflected in the PES. The robust link we have established with multi¬ 
pole response thus implies a relation between ionisation pattern and multipole 
response. This is certainly an interesting connection worth being investigated. 
Work along that line is in progress. 
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